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Nature of chiral spin liquids on the kagome lattice 
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We investigate the stability and the nature of the chiral spin liquids which were recently uncovered in extended 
Heisenberg models on the kagome lattice. Using a Gutzwiller projected wave function approach - i.e. a parton 
construction - we obtain large overlaps with ground states of these extended Heisenberg models. We further 
suggest that the appearance of the chiral spin liquid in the time-reversal invariant case is linked to a classical 
transition line between two magnetically ordered phases. 


Introduction — The quest for quantum spin liquids [1] 
is currently a very active endeavour in condensed matter 
physics. This elusive state of quantum matter comes in various 
forms and is theoretically intensely studied, however was dif¬ 
ficult to pin down in computational studies of realistic quan¬ 
tum spin Hamiltonians and hard to characterise unambigously 
in experiments on quantum magnets. 

The S' = 1/2 Heisenberg antiferromagnet on the kagome 
lattice has emerged as one of the paradigmatic systems where 
quantum spin liquid phases are expected. A plethora of the¬ 
oretical proposals have been put forward, ranging from va¬ 
lence bond crystals [2-6], algebraic spin liquids [7-10], Z 2 
spin liquids [11-17], to chiral spin liquids [2, 18, 19]. De¬ 
spite tremendous theoretical and computational progress [ 20 - 
32], the true nature of the ground state and the low-lying ex¬ 
cited states of the nearest neighbour Heisenberg model on the 
kagome lattice is still not settled completely. 

Chiral spin liquids (CSL) are a particular family of spin liq¬ 
uids in which time-reversal symmetry (TRS) and parity sym¬ 
metry are (spontaneously or explicitly) broken [33, 34]. The 
scalar chirality {Si ■ {Sj x Sk)) is non-zero and uniform and 
manifests the breaking of time-reversal and parity symmetries, 
analogous to the presence of an orbital magnetic field. In a fa¬ 
vorable situation the breaking of these symmetries could con¬ 
ceivably lead to a spin analogue of the Fractional Quantum 
Hall Effect, although other types of ground states are pos¬ 
sible as well [35, 36]. Historically Kalmeyer and Laughlin 
envisioned such a scenario by considering lattice versions of 
the bosonic z/ = 1/2 Laughlin wave function as candidate 
ground state wave functions for the triangular lattice Heisen¬ 
berg model [37, 38]. 

In two recent papers [39, 40], two forms of chiral spin liq¬ 
uids have been discovered, which are stabilised away from 
the nearest neighbour Heisenberg model upon adding further 
neighbour Heisenberg interactions or scalar chirality terms to 
the Hamiltonian. Both studies numerically demonstrate the 
required ground state degeneracy and characterize the under¬ 
lying topological order by computing the modular matrices. 

This breakthrough lays the foundation for further investiga¬ 
tions of chiral spin liquids. Several pressing, important ques¬ 
tions arise: i) are the two chiral spin liquids phases distinct 
or are they related ? ii) is there a simple physical (lattice- 
based) picture or a variational wave function that describes 
the chiral spin liquid ? iii) what is the ’’raison d’etre” of these 
chiral spin liquids, i.e. why are the chiral spin liquids stabi¬ 


lized for the two reported Hamiltonians ? Can we come up 
with some guiding principle which will allow to stabilise CSL 
on other lattices ? In the following we will address each of 
these questions. In short we find that the two chiral spin liq¬ 
uids are indeed connected. We then demonstrate that appro¬ 
priate Gutzwiller projected parton wave functions can have 
large overlaps with the numerically exact ground states of the 
studied microscopic models. And finally we show that one lo¬ 
cation of the chiral spin liquids in parameter space coincides 
largely with a transition line in the phase diagram of the cor¬ 
responding classical model. The classical transition line lies 
between coplanar q = 0 magnetic order and a chiral, non- 
coplanar magnetically ordered phase (cubocl [18]). 

Model — We will consider the following Hamiltonian 
which unifies the two models studied in Refs. [39, 40]: 

H=Ji^S^-Sj+J2 X S^-Sj + 

J3 ^ S,-S, + J^ Y. S,-{SjxSk). 

This model includes first, second and third nearest neigh¬ 
bour Heisenberg interactions with coupling constants Ji, J 2 , 
J 3 as sketched in Fig. 1 . The third nearest neighbour Heisen¬ 
berg interactions are only considered across the hexagons. 
While these interactions preserve TRS and all the discrete lat¬ 
tice symmetries of the kagome lattice, the additional three- 



FIG. 1. Sketch of the kagome lattice and of the different interac¬ 
tion terms of the Hamiltonian (1). Heisenberg interactions between 
first, second and third nearest neighbour are considered. The third 
nearest neighbour Heisenberg interactions are only considered across 
the hexagons. Three-spin scalar chirality interactions, breaking time- 
reversal and parity symmetries, are also considered on grey shaded 
triangles. 
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FIG. 2. Excitation Spectra from Exact Diagonalization. Different 
symbols and colors correspond to different momentum/pointgroup 
symmetry sectors. We use the cluster geometries and notation ex¬ 
plained in Ref. [27] (a) Effect of term on spectrum on the 30 sites 
cluster. The four-fold degeneracy of the ground state is lifted to a 
two-fold degereracy which corresponds to one sign of the scalar chi¬ 
rality. (b) Scan across the classical transition line for = 0 on 
the 36b sites cluster. The four-fold degereracy of the CSL is only 
present close to J 2 = J 3 (yellow shading), (c) Energy spectra for 
J 2 = J 3 = 0.4, = 0 and various system sizes Ns and geome¬ 

tries. Turquoise rectangle'. (0,0) [T] momentum, even under 180° 
rotation. Blue up triangle'. (0, tt) [M] momentum, odd under 180° 
rotation. Red down triangle: (0,0) [T] momentum, even under 180° 
rotation, odd under reflection. 


spin scalar chirality interactions on the triangles parametrized 
by break explicitly TRS and spatial parity. Note that 
Hamiltonian (1) features SU(2) invariance in spin space. For 
simplicity we will set Ji = 1 in the following. 

In Ref. [40] a CSL phase was found for O.OStt < 
arctan|(^| < tt/2 and J 2 = J 3 = 0. In this case, 
TRS is explicitly broken. Interestingly a two-fold degener¬ 
ate ground state was found, which furthermore exhibits the 
expected modular data and entanglement spectrum for a topo¬ 
logically ordered chiral 1 / = 1/2 Laughlin state-like phase. 
On the other hand in Ref [39] a chiral spin liquid with 
spontaneous TRS breaking was discovered for = 0 and 
0.2 < (J 2 = J 3 )/Ji ^ 0.7. Here the ground state degeneracy 
is four, which can be understood as arising from two copies of 
opposite chirality of a two-fold degenerate 1 / = 1/2 Laughlin 
state. Unlike several topological phases as Toric code [41] and 
double-semion [42] phases that also have a four-fold ground 
state degeneracy, we will show that in this case time-reversal 
symmetry is spontaneously broken. 

Energy spectroscopy — To investigate the persistence of 
this chiral spin liquid at the thermodynamical limit, we stud¬ 
ied the model for J 2 = J 3 = 0.4 and = 0 up to 42 
sites. The low-energy spectra for different system sizes are 
shown in Fig. 2 c). While the energy splitting between the 
four ground states has a non-monotonous behaviour, the en¬ 


ergy gap between the four lowest energy states and the hfth 
one increases with the system size. Moreover, the ratio of the 
energy splitting to the energy gap decreases with the system 
size, this tends to indicate that this phase is indeed realized 
at the thermodynamical limit. It is also important to notice 
that the momentum sectors involved in the four-fold degener¬ 
ate manifold depend on the cluster shape and can be predicted 
in complete analogy to the Fractional Quantum Hall and Frac¬ 
tional Chern insulator states [43, 44]. 

In Fig. 2 a) we investigate the energy splitting of the four 
ground states as we switch on a hnite coupling. At = 0 
the long-range order in the spin chirality is spectrally encoded 
in the presence of two states per topological sector, where the 
two states have to be at the same momentum, but differ in 
the spatial reflection quantum number (if the sample allows 
this symmetry). As is shown in Fig. 2 a), the two states per 
sector split very rapidly upon switching on 7 ^ 0. We can 
understand the action of regarding the scalar chirality in 
analogy to the effect of a longitudinal magnetic held on the 
two degenerate ground states in a ferromagnetic Ising model 
in the ordered phase, where the magnetic held immediately 
selects one of the two ordered states. As we show later based 
on overlaps, the chiral spin liquid thus selected by is of the 
same type as the one stabilised in the Ji — model alone, and 
is connected to the TRS symmetric situation in the absence of 

Finally in Fig. 2 b) we investigate the effect of a deviation 
from the J 2 = J 3 condition (in the absence of J^) by Ax¬ 
ing J 2 = 0.5 and varying J 3 . One observes that the four¬ 
fold ground state degeneracy is rapidly lifted when J 3 devi¬ 
ates more than about 0.05 ~ 0.1 from 0.5. Interestingly the 
line 0 < J 2 = >73 < 1 is the classical transition line be¬ 
tween a magnetically ordered q = 0 ground state for J 3 < J 2 
and the non-coplanar magnetically ordered cubocl phase for 
•h > J 2 [45]. Below we will show that also the overlaps 
with the variational wave functions are large only in the direct 
vicinity of this classical transition line. A deeper understand¬ 
ing of the classical ground state conflgurations on that line and 
of the effect of quantum fluctuations on that manifold might 
thus lead to an identification of the crucial ingredients required 
to predict and uncover chiral spin liquids in TRS Hamiltoni¬ 
ans on different lattices. We note in passing that the explicitly 
TRS breaking Hamiltonian (1) with =/= 0 can be consid¬ 
ered as a truncated version of a parent Hamiltonian for the 
CSL constructed in Ref. [46], similar to the spin Hamiltonian 
on the square lattice considered in Ref. [47]. 

Parton construction and overlaps — As stated earlier 
on, the CSL can be considered as lattice analogues of the 
bosonic v = 1/2 Laughlin state. In recent years substan¬ 
tial activity focused on realizing such states on fractionally 
filled Chern insulators, so called Fractional Chern Insulators 
(FCI) [43, 48, 49]. It is thus a natural question whether the 
CSLs under consideration might also have such an interpreta¬ 
tion. The natural bosonic v = 1/2 FCI state on the kagome 
lattice [50, 51] however does not have the correct magnetisa¬ 
tion since it corresponds to magnetisation m/wsat = ±2/3 
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FIG. 3. Band structure of the [7r/2,0] - model used to generate 
the pair of model states with positive scalar chirality, a) geometry 
of the unit cell. The mean-field parameters Xij are chosen such that 
there are 7r/2 fluxes through the triangles and no flux through the 
hexagons, b) Brillouin zone of the model (shaded) and conventional 
first Brillouin zone of the kagome lattice (hexagon), c) Band struc¬ 
ture of the model along the path between high symmetry points in 
the Brillouin zone as drawn in b). The bands are separated by a finite 
gap. Each solid (dashed) band carries Chem number —1 (+1). 


instead of the required m = 0 [52]. 

In the absence of a simple FCI candidate wave function we 
pursue an alternative approach, based on a parton construc¬ 
tion. In order to understand and classify the different spin 
liquids a generalized construction scheme called parton con¬ 
struction has been introduced by Refs. [53-58] - see [59] for 
an introduction. The main idea of this technique is to split 
up each spin operator Si at site i into two fermionic parton 
operators according to 

Iff 

Note that by introducing these operators the Hilbert space is 
enlarged due to the possibility of doubly occupied or vacant 
sites. Substituting the parton operators for the spin operators 
and performing a mean-held approximation by introducing 
mean-field parameters Xij — Tha yields (ignoring 

constants) a tight-binding model of type 

f?mean = ^ {xijclrjCja + h.C.^ . (3) 

Several of these models have been investigated for the 
kagome lattice [2, 7-10]. Here we focus on nearest neigh¬ 
bour Xij only and the norm is chosen to be \xij\ — 1- 
Physically different states can be created by choosing Xij 
such that different magnetic fluxes thread the triangles and 
the hexagons of the kagome lattice. Amongst these states 
we consider states whose parent mean-field models have uni¬ 
form ±7r/2 flux through the triangles and zero flux through 


the hexagons [2, 7, 8]. To do so a magnetic six sites unit 
cell is needed instead of the standard three sites unit cell of 
the kagome lattice. In the following we will call these the 
[±7r/2,0] - models. 

The unit cell geometry, Brillouin zone and band structure 
of the [7r/2,0] - model are shown in Fig. 3. All six bands 
have non-zero Chern numbers as indicated in Fig. 3c). To 
obtain a = 0 model state, the three lowest bands are com¬ 
pletely filled both for spins up and spins down and an exact 
Gutzwiller projection is applied to project onto the physical 
spin subspace. As the filled bands are separated by a finite 
gap from the empty ones, the spin-spin correlations after pro¬ 
jection are expected to decay exponentially with distance, and 
thus describe a spin disordered state. The Chern number of 
the filled bands for the [7r/2,0] - model ([—7r/2, 0] - model) is 
— 1 (+1). The [7r/2,0] - model ([—7r/2, 0] - model) yields a 
positive (negative) scalar chirality expectation value for every 
basic triangle. 

On the torus there are two independent non contractible 
loops. Some of the gauge choices which leave the flux through 
the triangles and hexagons invariant, correspond to different 
fluxes through these torus loops. Threading flux through these 
loops corresponds to a Laughlin flux insertion. Thereby dif¬ 
ferent topological states can be generated. These states cannot 
be distinguished by local observables and therefore are de¬ 
generate for local Hamiltonians in the thermodynamic limit. 
For the chiral spin liquid a two-fold topological ground state 
degeneracy is expected. Thus, by threading different fluxes 
through the torus we should only be able to create a two- 
dimensional space. We numerically computed the Gutzwiller 
projected wavefunctions (GPWFs) of the [±7r/2, 0] - models 
with a fixed gauge. In order to construct the topological part¬ 
ners of these states we additionally thread fluxes through the 
torus as explained in the previous section. We checked that for 
each of the [±7r/2, 0] - models, only two linearly independent 
states can be only constructed as expected for a CSL within a 
numerical accuracy of 10“^, similar as in Ref. [60]. 

We compare now these four model states with the ground 
state |'0ed) of the Hamiltonian (1) obtained using Exact Diag- 
onalization. We choose the overlap of the ground state 
wave function with the four model states as our figure of merit; 


KV'edIV': 


GW! 


(4) 


Overlaps of the GPWFs with the ground state of the Hamilto¬ 
nian (1) for different parameters on a = 30 sites sample 
are shown in Fig. 4. The overlaps of our model state with the 
ground state wave functions of the model of Ref. [40] where 
J 2 — J 3 = 0, Ji = cos 9 and = sin0 are shown in 
Fig. 4(a). We found that overlaps for between 0.1 and 1 
range from 0.62 to 0.97. The overlap of the two GPWFs of the 
[—7r/2,0]-model are by orders of magnitude larger than those 
of the [+7r/2,0]-model. For between —1 and —0.1 the 
overlaps are exactly the same within numerical precision as 
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FIG. 4. Overlaps Oqw of GPWFs with ground states from Exact Diagonalization for J 2 = J 3 = 0 and Ji = .y /1 — on a 30 sites lattice. 
The overlaps of the GPWFs of the ± 7 r/ 2 -models are symmetric under changing the sign of J^. The maximum overlap is equal to 0.97 and is 
reached for = ±0.7. The CSL phase extends almost up to the Heisenberg point. 


for between 0.1 and 1 but the role of the GPWFs from the 
[± 7 r/ 2 ,0]-model and [— 7 r/ 2 , 0]-model are exchanged. This 
is expected since the model with negative should have a 
positive scalar chirality and therefore only little overlap with 
the variational states form the [— 7 r/ 2 , 0 ]-model with negative 
chirality and vice versa. 

For the time-reversal symmetric model with = 0, our 
variational wave functions have substantial overlap only close 
to the line J 2 = J 3 , in agreement with the energy spec¬ 
troscopy results discussed above [Fig. 4(b)]. In this region 
the overlaps reach up to 0.72 for Ns = 30. 

As can be seen in Fig. 4 c) for — 0.3 and for = 
0.6 (not shown), the region of the CSL broadens signihcantly 
when is increased from zero. For = 0.3 (resp. = 
0.6) the overlaps on the classical transition line for J 2 = J 3 
between 0 and 0.7 range from 0.8 to 0.9 (resp. from 0.85 to 
0.95). 

Conclusion — We showed that the two recently found re¬ 
alizations of chiral spin liquids on the kagome lattice [39, 40] 
are indeed related and can be described by Gutzwiller pro¬ 
jected parton wave functions. This yields an intuitive micro¬ 
scopic picture of the CSL phase stabilized in these models. 
The ansatz wave functions we chose have been shown to de¬ 
scribe a CSL on the kagome lattice [2, 7, 8 ]. We constructed 
a pair of Gutzwiller projected parton CSL wave functions for 
each sign of the scalar chirality. We suggested that these states 
describe the CSL ground state found on the kagome lattice. To 
prove that indeed these wavefunctions describe the novel CSL 
phases found in Refs. [39, 40] we computed overlaps of these 
variational wave functions with the ground state wave func¬ 
tions computed by Exact Diagonalization. Substantial over¬ 
laps were found in regions of the phase diagram where the 
CSL is expected. By further investigation of excitation spec¬ 
tra, we showed that the CSL phase in Ref. [39] is only present 
on the transition line between a chiral cubocl and a coplanar 
q = 0 phase of the classical phase diagram [45]. This could 
serve as a guiding principle for hnding CSL phases in other 


models and on others lattices. Being related to the Laughlin 
state, these states should exhibit anyonic excitations. Their 
investigations will be pursue in a future work. 

Note added — While completing the present manuscript 
we became aware of parallel work reaching similar conclu¬ 
sions using complementary methods [61, 62]. 
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